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Abstract 

We derive from the Navier-Stokes equation an exact equation satisfied by the dissi- 
pation rate correlation function, < e(x + r, t + r)e(x, t) >, which we study in the equal 
time limit. We exploit its mathematical similarity to the corresponding equation derived 
from the 1-dimensional stochastic Burgers equation to show that the main intermittency 
exponents are \x\ = 2 — and \xi — 2z^ — £4, where the £'s are exponents of velocity 
structure functions and £4 is a dynamical exponent characterizing the 4th order structure 
function. We discuss the role of sweeping and Galilean invariance in determining the 
intermittency exponents. 

PACS numbers: 0.5.45.+b, 47.10.+g 

1 Introduction 

The correlations of the local energy dissipation rate per unit mass, e(x, t), whose behavior in 
the inertial range can be written as 

< e(x + r)e{x) >~ < e > 2 (r/L)~^ . (1) 

play an important role in the study of high Reynolds number turbulence jl], H|. The reason 
is that the algebraic decay of the correlations characterized by a positive value of [i signals 
intermittent behavior. In Eqn. (1) < e > is the mean value of the dissipation rate per 
unit mass, r belongs to the inertial range, and L is the integral scale or system size. We 
have assumed that the turbulence is homogeneous and isotropic. When /j, = there are no 
fluctuations at large distances and scaling holds. A positive intermittency exponent /i signals 
the breakdown of Kolmogorov scaling and, experimentally this is indeed the case; there now 
appears to be agreement [[| on a value of /i = 0.25 ± 0.05. 

The standard description of the statistical properties of turbulence in terms of velocities 
leads to the introduction of (longitudinal) structure functions, S p (r), whose behavior in the 
inertial range is given by 

S P {r) = < [(u(r) - u(0)] • r] p > ~ r & . 



A scaling argument (of which a more sophisticated version is the refined similarity hypothe- 
sis of Kolmogorov and ObukhovQ) has been employed to relate the exponent fi character- 
izing dissipation to the exponent, which describes the behavior of the 6th order struc- 
ture function. In the naive scaling argument, since e has dimensions of velocity 3 /length, 
one expects that e, when it occurs in correlation functions, behaves as (Au) 3 /r, where 
Au(r) = [u(x + r ) — u(x)] ■ r is the longitudinal velocity difference. Therefore, the dissi- 
pation correlation in (1) can be expected to decay as SQ(r)/r 2 . This reasoning leads to the 
well-known identification 

M = 2-Ce • (2) 

Experimental measurements yield ^6 ~ 1.8 (its scaling value is 2); the value of n which is then 
obtained from Eqn. (2) is consistent with the experimental result of u = 0.25 ± 0.05 quoted 
above. Because of this agreement many discussions concerning the intermittency exponent 
have been limited to a discussion of the non-scaling behavior of although other expressions 
for a have been proposed. In particular, the relation 

H = 2C 2 - C4 (3) 

has been proposed. |||, ||] It leads to a value for a of // ~ 0.1 using data in the literature. || 
For r « L , the larger value of u would dominate. 

Indeed one should generally expect that in an appropriate equation satisfied by the dis- 
sipation correlation different terms will lead to different possible values for a from which the 
dominant behavior in the inertial range can be extracted. As far as we know a systematic 
treatment of this question does not exist. The aim of this paper is to provide one. We present 
an exact equation satisfied by the dissipation correlations, analyze it, and in particular, show 
how the two aforementioned values of fi arise. We should stress here that the only way we 
have found to derive the equations for the (experimentally measured) equal time, spatial 
correlations is to start from correlations in both space and time, and then impose the equal 
time limit. 0] 

We will treat both the Navier-Stokes equation and the one-dimensional stochastic Burgers 
equation. || The mathematical structure of equations derived from the latter, while simpler 
algebraically, is similar to that encountered for the Navier-Stokes equation, as we have pointed 



out previously. ||, |TT], [TJ| The stochastic Burgers equation thus provides a gateway to 
understanding issues of fully developed turbulence and our results confirm the usefulness of 
this approach. It enables one to go beyond generalized scaling arguments which essentially 
embody the correctness (up to intermittent effects) of the original Kolmogorov argument. 

In Section 2 of the paper we present the exact equations for dissipation correlations, 
both for the Navier-Stokes equation and the stochastic Burgers equation. The details of the 
straightforward derivation are relegated to Appendix A. Our analysis discussed in Sec. 3 
shows that there are two main contributions in the inertial range: the dominant one charac- 
terized by the exponent a = 2 — (q, and a second which depends on the dynamical behavior of 
the 4th order structure function. The mathematical similarity of the analysis for the Burgers 
problem and that for the Navier-Stokes equation is discussed. In section we turn to a 
discussion of some technical points, in particular, the role played by Galilean invariance of 
the two fundamental equations, and its relationship to "sweeping" effects. A summary of our 
results and conclusions follow. Technical points are taken up in Appendices B, C and D. 
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2 The fundamental equations and their derivation 



In addition to the Navier-Stokes equation we will also investigate, as mentioned in the Intro- 
duction, the one-dimensional, stochastic Burgers equation, || describing a compressible fluid, 
given by 

du/dt + udu/dx = vd 2 u/dx 2 + rj(x,t) (4) 

where v is the kinematic viscosity, and ry a long-ranged, Gaussian random noise with zero 
mean and correlations in k space given by 

< fj(k, t)fj{k\ t') = 2D(k)5 k+k , fi 6(t - t') (5) 

where D(k) = Dol^l^- The exponent (5 determines the range of the forcing and we focus on 
j3 € (0, —1]. We have shown previously || 10 1 that for these values of (5 the system exhibits a 
rich multifractal structure, going from a near scaling regime for (5 close to zero to an extreme 
multicritical regime characterized by the presence of strong shocks, with Q n — > 1 for n > 3 
when (5 — > —1. 

We will derive from the stochastic Burgers equation and analogously in the Navier-Stokes 
case, an equation satisfied by the correlations of the local dissipation rate. We emphasize 
that it is convenient to derive an equation for the time-dependent quantity 

C(r,r) =< e(x + r,t + r)e(x,t) > (6) 

where the local dissipation rate is defined by 

e = ^(diUj + djUi)(diUj + djUi) (7) 

in the Navier-Stokes case and by e = v{du/dx) 2 in the one-dimensional Burgers problem. We 
have employed the notation di = d/dxi. The equal time correlation function can be obtained 
by taking the limit r —* 0. 

Denote the space-time coordinates of two points by x and x' = x+r and t and t' = t+r and 
the various fields by primed and unprimed variables, for example e' = e(x' , t'), e = e(x, t) and 
similarly for the velocities and other fields. We consider the Navier-Stokes equation driven 
by a random force fj(k, t) concentrated around a small wavevector; this enables us to define 
average values most simply. It is straightforward to derive (see Appendix A) the following 
equation describing the spatio-temporal behavior of the correlations of the dissipation rate: 



< e'e > = —d T < eu' ■ u' — e'u ■ u > + -cV < eu' ■ u'u'a — e'u ■ uiij > 

A A ■> J 



—d T < eu ■ u' — e u ■ u > + -< 
4 4 

+— V r < e u ■ u + eu' ■ u' > +— V r V r . < e muj + eUiUj > 

+ - < e'u ■ fj + eu 1 ■ 7]' > < ^V u i ~ e P >u i > (8) 

where p is the pressure divided by the mass density. This exact equation is one of the fun- 
damental equations in this paper. 

Similarly one can derive (see Appendix A) the corresponding equation for e(x, t) in the 
Id stochastic Burgers equation. We will denote the two space-time points by x\ = x + ^r, 
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t\ = t+ gT, X2 = x— 2?", and £2 = t— 2 t respectively and use the notation u\ = u{x\,t\), etc., 
for the various fields. The equation satisfied by C(r,r) =< eie 2 > is =< e(xi, ti)e(x 2 , t 2 ) > 

< ei£2 > = t^t < ei^2 2 - £2^i 2 > +7r<9 r < ei^2 3 - e2""i 3 > 
4 b 

^ o2 2 2 

+-d r < ei« 2 + £2^1 > 

+7: < £iU2V2 + ^uiVi > ■ ( 9 ) 



We draw attention to the close mathematical similarity of the two equations, in particu- 
lar, to the first two terms in each equation, schematically d < eu 2 > jdr and d < eu 3 > /dr, 
which yield the dominant contributions to the behavior of the dissipation correlation. The 
second line in each equation contains terms proportional to v and these do not contribute 
in the large Reynolds number limit. We emphasize that the z^— dependent terms that are 
multiplied by singular operators have to be handled carefully to extract the innocuous terms 
in the second line (see Appendix A). The difference in the last lines is evidently due to the 
pressure term in the Navier-Stokes equation. We hasten to add that physically the nature of 
turbulence is not the same in the two systems. For example, vortex stretching is important 
in the three-dimensional problem. In addition, the energy flux through a wavenumber ko, 
n(fco), is proportional to (Dq/ L)k^ + ^ in the stochastic Burgers problem while it is a con- 
stant in the Navier-Stokes case. Nevertheless, the mathematical similarity permits a parallel 
analytical investigation. Thus the main features of Eqn. (g), namely the expressions for the 
dominant intermittency exponents quoted in the Introduction, can be extracted by learning 
to analyze the equation in the one-dimensional stochastic Burgers equation and employing 
a similar strategy in the Navier-Stokes case. These equations for the dissipation rate cor- 
relations depend crucially on including time derivatives and the variation of the temporal 
fluctuations contribute to the potential intermittent behavior of < ee' >. 



3 Analysis of the fundamental equation in the Id Burgers case 

We will begin with the mathematically simpler one-dimensional case; in the rest of the paper 
we restrict our attention to the limit r — > 0. In this limit the last term in Eqn. (||) which 
depends on the random force is equal to the product < e\ >< £2 >, apart from terms that 
vanish in the limit v — > 0. This is shown most easily in k space, given the correlation of the 
random force in Eqn. (||), using Novikov's result for Gaussian random processes[13| and the 
result < e >= jp J2k D(k). This leads in the r — > limit to 

C(r) =< eie 2 > - < e > 2 = V -[ d 2 r < e 2 V + em 2 2 > +d r S 2 (];d 2 S 3 - vd?.S 2 ) } 

4 D 

+\d r < e 2 «i 3 - eiu 2 3 > +\d T < e 2 u 2 - e x u 2 2 > .(10) 
o 4 

All the terms on the right-hand side are evaluated in the r — ► limit. 
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3.1 Outline of the argument and main results 

We analyze the fundamental equation in the Burgers case, Eqn. (|lOl), for r in the inertial 
range and extract the dominant terms. In the zero-viscosity limit it is easy to see that the 
first line in Eqn. ( |1C)| ) does not contribute: clearly, terms that do not depend explicitly on 
e are finite and hence, are negligible as v — > 0. Note that in contrast to the case of fully 
developed turbulence where < e > is a constant, in the stochastic Burgers equation one has 
< e >oc 5~ 1-/3 where 5 is an inner cutoff scale (shock thickness) and since v oc 5 l ~^^ the 
other terms proportional to v vanish for j3 < 0, fl4}| 

The crux of the argument for extracting the possible values of ji depends on the obser- 
vations that the dominant behavior of < eiU2 3 — e2^i 3 > is determined by dSe(r,r)/dr = 
d/dr < ui — U2) 6 > and that of < ei«2 2 — ^2U\ 2 > primarily by dS^/dr (we discuss this 
later) and thus the dominant singular behavior of the dissipation correlation contains 

d 2 S 6 J fd 2 SA 
and 



dr 2 \ <9t 2 / 

\ / r=0 

Since Se(r) ~ r^ 6 in the inertial range the first term immediately yields the well-known 
intermittency exponent 

Ml = 2 - Ce (11) 

for the large distance behavior of C(r). 

The second term leads to an intermittency exponent 

fi 2 = 2z~ 4 - C 4 (12) 

where the exponent 2z^ determines the dynamical behavior of the second derivative of the 
fourth order structure function when r — > 0. If one assumes naive dynamic scaling |l5[| 



one has £4 = Z4 and thus 



S A (r, T ) « \r\^ S 4 (\r\/\rn (13) 



lii = 2z 4 - U ■ (14) 



We will show later that similar results obtain in the Navier-Stokes case also. Accepting this 
for now and using the fact that in the scaling limit, there is only one dynamical exponent 
z = Z4 which has the same numerical value as £2 = 2/3 in the Navier-Stokes case[jl]], the 
second intermittency exponent can be written as 

[12 = 2( 2 - C4 

which is Eqn. (^). This is precisely the result obtained in reference § using scaling ar- 
guments for static expressions. Our analysis shows that the expression in Eqn. (||) is an 
approximation, in the scaling limit, of what is actually a dynamical result. In the Navier- 
Stokes case if Kolmogorov scaling holds fix = fJ-2 = 0. Returning to the stochastic Burgers 
equation, for small (3 where all the relevant structure functions scale fJ-i = fJ-2 = 2 + 2/3; in the 
limit (3 —* — 1 where Q2 = 2/3, the Kolmogorov value in 3d turbulence, the scaling values of 
both [i\ and fi2 vanish as in the Kolmogorov case. Note, however, that n\ > 1 in the range 
P E (0, 1]. As (3 — > — 1 and multiscaling occurs for p > 4 , \i\ = 2 — £g approaches unity 
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from above and dominates the behavior of C (r) . [fig] We draw attention to the fact that for 
the dynamical exponent in Eqn. (12) we have not used the value of the kinematic exponent 
z = 1 that arises from sweeping but the intrinsic dynamical exponent. In Sec. 3.4 we provide 
a more detailed justification for this. 



3.2 Contribution of the < exw| — e 2 u\ > term 

We will now argue that < eiu\ — e 2 uf > behaves as d r Se for r = as claimed above. In 



an earlier paper 1 11] we have presented the equations satisfied by the (low-order) equal-time 



structure functions; in particular, the sixth-order structure function, Sq, satisfies the equation 

d g6( X = 0) = 3l/ ^ 2 S 5 (r) - 30 < (ei + e 2 )(«i - u 2 f > + y < («i - u 2 )\ m - m ) > 

= 3ud 2 S b {r) - 30[< (ei + e 2 )(«i - u 2 f > - < e x + e 2 > 5 3 ] + 5S 3 (r)dS 3 /dr 

(15) 

In obtaining the second equality, the noise term has been simplified using Novikov's theorem 
and the von-Karman-Howarth relation. The contribution S3 dS 3 /dr that arises from the noise 
term can clearly be identified as the scaling contribution to Sq which yields for the exponent 
describing the inertial-range behavior £g = ^£3 = —2/3. The v— dependent contribution 
vanishes as v — > since it is multiplied by a term which is finite in the inertial range and thus 
we have, generically, the result that terms of the form < e\u\ > are proportional to dS%/dr. 
This immediately yields the result m = 2 — £g as stated earlier. We emphasize the feature 
that the only term that can possibly lead to non-scaling behavior of S§ is the term which is 
a product of the dissipation rate and velocities. This result, derived from the exact equation 
for Sq, agrees with the field-theoretic point of view in which the non-scaling behavior arises 
from short-distance singularities which are determined by the limit v — > and the existence 
of an energy cascade. 



3.3 Contribution of the < €iu\ — e 2 u\ > term 

In contrast to the contribution of d T < eiu 2 — e 2 u\ > to the dissipation correlation function, 
the term d T < ei« 2 — e 2 u\ > | T= o is an intrinsically dynamical object. In order to study it 
we consider the temporal evolution of S^r, r) =< (u(xi, t±) — u(x 2 , t 2 )) 4 > where r = t\—t 2 
and r = x\ — x 2 . It can be shown to satisfy the equation 

dSi/dr = - ^d r < (tta -iii) 4 (n 2 + m) > + 2 < (m - u 2 ) 3 (r?i + m) > 

-6 < (ei - e 2 ){u 2 - Ul ) 2 > . (16) 

The derivation of the general equation satisfied by the time evolution of the generating 
function for S q (r,r) is given in Appendix B. The above equation follows directly from it as 
a special case. This equation contains a term of the generic form < eu 2 > whose derivative 
with respect to r occurs in the basic equation for < e\e 2 >. The analysis of this equation is 
somewhat more subtle; however, once again we expect that the non-scaling behavior of the 
left-hand side is determined by correlations that involve products of the dissipation rate and 
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velocities, i.e., the last term on the right-hand side. As in the previous subsection we expect 
the noise term to yield a scaling result: to see this we first note that in the limit r — > the 
noise term yields a discontinuity since it gives a different contribution depending on whether 
r — > + or r — > 0~ . A careful analysis yields, for r > 0, a contribution 

which dominates the behavior of d T S^\ T= Q. In the range of (3 we are interested in both S2 



and S3 exhibit scaling behavior. [17] Thus, the noise term again yields scaling behavior for 



the dynamics of the first derivative in the r — > limit. The first term in Eqn. (16) ensures 
Galilean invariance of the equation and contains only sweeping effects as will be discussed 
below. Therefore, the behavior of terms of the form \d T e\u^\ T= Q in the equation for < eie2 > 
will reflect the non-trivial singularity in d^S^ | r =o> which then yields the second intermittency 
exponent, expression (12), namely \xi = 2z^ — (4. 



3.4 Technical remarks 

In this section we discuss some of the technical points that were not addressed in the ear- 
lier sections. There are two related difficulties: the first is the more obvious one; while our 
equation for < e\€2 > (see Eqn. (JToj) ) involves velocities, our equations for d r SQ and 8-,-Sa 
involve velocity differences and the precise combinations of the terms do not match exactly. 
Second, if we carry out a straightforward analysis of the term d T < eiii| — £2^2 > we obtain 
terms of the form d% < (u\ — "U2) 4 (^i + ^2) >• A naive factorization of the expectation 
value leads to an exponent 2 — ^4 and not the value of fi2 claimed earlier. In particular, 
this corresponds to the kinematic value Z4 = 1 in the expression for [12 (Eqn. flli])) and not 
to the dynamical value of 24, which in the scaling limit is numerically equal to Qi. While 
we have not been able to provide a mathematically complete analysis of all these points we 
present below convincing arguments in support of our results based on general considerations. 

The first observation we make is that while < £162 > is manifestly Galilean invariant 
the invariance of the right hand side of Eqn. (|f^) is not evident; it can be checked explicitly 
(see Appendix C). Therefore, we expect a (Galilean-invariant) operator product expansion 
analysis of the terms of the form d T < exu| — £2^1 > to pick out precisely the terms that 
we have used in our analysis since the terms will involve Galilean invariant combinations of 
operators such as, for example, velocity differences. Physically, the potentially troublesome 
terms, such as 2 — (4 mentioned above, arise from the effects of sweeping. For example, we 
have discussed in a previous paper (T^j how behavior such as 

< (ui - u 2 ) A {u\ + uf) >~< («i - u 2 ) 4 >< (u\ + uf) > 

arises. For nonzero values of r such a Galilean non-invariant term is known to contain ef- 
fects of sweeping. That our formalism gives rise to such effects is not surprising since we 
have adopted an Eulerian point of view. As pointed out by Tennekes|l8|] large-scale energy 
containing structures advect inertial-range information past an Eulerian observer. This al- 
lows one to note that the dominant term in < {u\ — U2) 4 (n| + n|) > is proportional to 
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< u 2 >< (u\ — U2) 4 > where we have exploited the independence of the large-scale flow 
characteristics encapsulated in < u 2 > and the information about the inertial range con- 
tained in < u\ — U2) 4 >■ Good numerical support for the correctness of such an argument 



has been presented in ref. [12]. In other words, in the Eulerian picture the dynamical be- 
havior of the structure functions at finite values of the time difference r is dominated by 
sweeping effects which are characterized by an effective dynamical exponent of z = 1. (The 
second intermittency exponent, fi2, shown earlier to be 2z^ — Q would then reduce precisely 
to 2 — £4 as discussed in the previous paragraph. ) However here we are dealing with the limit 
t — > 0. Physically, equal time correlations do not show sweeping effects. Therefore, only the 
intrinsic, i.e., non-sweeping, dynamical behavior contributes to < eie 2 >■ In support of this, 



we also point out that we have shown previously |12j that intrinsic dynamical fluctuations 
determine the behavior of d T S2\ T =o an d yield the exponent Z2 = 1 + /3/3, while the correlation 
S^ir = 0,r), where one observes fluctuations at a given point as a function of time picks up 
the sweeping contribution of exponent z = 1. 

Experimentally, given the validity of Taylor's frozen turbulence hypothesis, especially in sit- 
uations with a large, externally imposed flow, one can determine equal time correlations by 
making one point measurements as a function of time. In fact, the transition from spatial to 
temporal variables depends on the kinematic exponent z = 1. Once, the equal time dissipa- 
tion rate correlations are thus determined, albeit approximately, their behavior is given by 
our analysis. 

One other term in C(r) not considered thus far is the r — » limit of 

^- < (ui - n 2 ) 3 (r/i + 772) > 
or 

. This term results from taking the djdr of Eqn. (|l6|). The term < (u\ — 112^ + f/2) > 
can be shown to be related to a fourth-order response function and a naive argument yields 
inertial-range behavior characterized by an exponent that is the scaling value of \%2- 

Thus, the behavior of C(r) in the r — > limit is determined by appropriate Galilean 
invariant combinations of spatial and temporal derivatives of structure functions, such as 
0t or 8 2 < ( Ul -u 2 ) 4 > +d T d r < (m -n 2 ) 4 (m +u 2 ) > +(l/4)<9 r 2 < (ux -u 2 )\u l +u 2 ) 2 >, 
and they give rise to the dominant exponents discussed above. Of course, we have not 
provided a first-principles determination of the dynamical exponent Z4, or for that matter, 
of ^6- We aimed only at establishing relations for the possible value of the intermittency 
exponent /i. 



4 Analysis of the Fundamental equation for the Navier-Stokes 
case 

We will follow in the three-dimensional Navier-Stokes case the strategy outlined earlier for 
the stochastic Burgers equation. The calculations are much more involved to carry out in 
the same detail. We will point out the mathematical correspondences between terms in the 
Navier-Stokes and the Burgers cases: we obtain the vectorial or tensorial generalizations of 
the scalar terms in the Id Burgers problem that are allowed by rotational invariance. For 
example, the analysis of the crucial < e(5u) 3 > term in Eqn. (|To| ) that gave rise to \x = 2 — (& 
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was based on the equation for d r SQ. Thus we are led to consider the corresponding equation 
for an appropriate sixth-order structure function in the Navier-Stokes case. We use the 
notation 



6u = u(x', t') — u(x, t) and Sff = ff(x', t') — fj(x, t) . 
for velocity and noise differences. We have derived the following equation 

— — < (8u ■ 5u) 2 5uj5uk > = ttV^ < [5u ■ 5u] 2 5u k > 
orj 2 

—4 < 5u ■ 5u 5u k [e + e'] > —6 < 5u ■ 5u 8iii[iki + e' ki ] > — 2 < 5ui5uj 5u k [£ij + e'ij] > 
< [Su ■ 5u] 2 5r]k + 45ukSu ■ 5u5u ■ Sff > 

_ 2 < [5u ■ 5uf d{p + p,) >-4<5u- 5u5 Uk 5 Ul d( - P + pr) > (17) 
or k dri 

In the preceding we have used the definitions for a (symmetric) tensor with the dimensions 
of the energy dissipation rate, 

ejfc = vd%Uj diU k . (18) 

The analogy of Eqn. ([l7| ) with the corresponding equation in the Burgers case (see eqn. 
(|l5|)) is striking. The first terms in both equations have the same form and vanish in the 
v — ► limit. The driving term < (5u) 4 5r] > in the Burgers problem generalizes to the two 
vector terms in the third line of Eqn. (|i~7|). The scalar term < (Su) 3 (e + e') >, which is crucial 
for identifying fj,\ in the Burgers problem, generalizes in the Navier-Stokes case to the three 
terms in the second line of Eqn. (|l7]) and involves the tensor generalization of e defined above. 
It is reasonable to expect that terms on the right-hand side generically of the form eu 3 which 
involve the dissipation rate lead to the intermittent behavior of J^- < [5u ■ 5u] 2 5uj5uk > 
since the dissipation rate contains the singular behavior of operators in the v — > limit. The 
other difference with the Burgers case arises from the pressure terms in the last line of Eqn. 
(|i~7|) and we do not expect it to yield the dominant intermittent behavior of the right-hand 
side since it does not depend explicitly on e; we expect it to yield scaling contributions. 



Straightforward tensor analysis shows that the left hand side of Eqn. (17) effectively includes 
the contribution of the sixth-order longitudinal structure function. We do not enter into the 
issue of whether longitudinal and transverse structure functions yield the same exponents; 
if they do differ as has been claimed |l^] then our results should be modified appropriately. 
The situation is simpler if as is argued by L'vov et q/. pO| there is a single exponent. 

We have also derived the equation for the temporal derivative of an appropriate fourth- 
order structure function. The details of the derivation are presented in Appendix D; the 
equation is given by 

d T < [Su ■ 5u} 2 > = - -d rj < (uj + u'^Su ■ 5u} 2 > +2([6u ■ 5u] 2 5u ■ (ff+ rf)) 

— 2 < (e — e')Su ■ Su > — 4 < (e ik — e' ik )5uiduk > — 
+2v{5u- 5uV 2 5p) 

- 4( 5u-5u^-(5p5ui) ) (19) 
on 
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Again we note the similarity to the structure of the equation in the Id case, Eqn. (U), 
which leads to the second intermittency exponent /U2- The first term on the right-hand side 
which ensures the Galilean invariance of the equation is the vector version of the first term 
on the right-hand side of Eqn. (|l~6l ) as is the noise term. The term that has the generic 
form < 5e(5u) 2 > which is a scalar in the Burgers case involves the tensor generalization of 
e defined earlier. It is easy to see that the term proportional to v is negligible in the high 
Reynolds number limit. The last term in Eqn. (|l{^) which depends on the pressure has no 
counterpart in the Burgers case; as argued earlier we do not expect this term to dominate the 
intermittent behavior of the left-hand side. On physical grounds it is reasonable to expect 
that the leading intermittent behavior of d T < [5u ■ 5u] 2 > is determined by terms that 
involve the dissipation rate on the right hand side, of the form < 5e(Su) 2 >, as in the Burgers 
case; thus the Galilean-invariant contribution of d T < eu' ■ u' — e'u ■ u > to the dissipation 
correlations in Eqn. (|8|) will be given by d 2 S4/ '3t 2 \ t= q. 

We see that the mathematical similarity of the various equations in the Burgers and 
Navier-Stokes cases allows us to carry out a similar analysis and deduce the same relations 
for the leading intermittency exponents in terms of the exponents that characterize the static 
and dynamic behavior of the velocity structure functions. We have discussed the implications 



of this already in Sec. 3.1 



5 Discussion of Results 

It is remarkable that the exponents which characterize the power law decay of the (equal 
time) dissipation rate correlations obey the same relations for the Burgers equation as for 
the Navier-Stokes equation. The leading exponent hi = 2 — (^q derives from the presence of 
8 2 Sq in the equation for dissipation fluctuations while the second exponent ^2 arises from 
the term d 2 S/i\ T= Q. Clearly, the origin of the second exponent (expression (14)) is dynamical, 
and coincides, in the Navier-Stokes case, with the static expression in Eqn. (||) (see Ref. [Q]) 
in the scaling limit, where the only dynamical exponent z is numerically equal to £2- 

We have derived the exact equation satisfied by C (r , r) and analyzed it to obtain our 



results. Recently, L'vov and Procaccia [21, 22] have found general fusion rules for multipoint 



correlations in the Navier-Stokes problem and used them to show that /x = 2 — £g; they also 
discuss the possibility of another scenario (see Eqn. (15) in Ref. |^lj) in which the static 
version of ^2, i-e, [12 = 2(2 — C4 can occur. Q Our approach finds these as possible terms from 
an exact equation and provides a dynamical expression for the second exponent. 

In our analysis we have had to disentangle sweeping effects from the dynamics of internal 
evolution, both of which are present in our Eulerian equations. We did this by exploiting 
the Galilean invariance of the internal dynamics, noting that sweeping which leads to a 
kinematic value of z = 1 yields effects that break Galilean invariance. It is clear that for the 
measurable quantity we are considering, namely dissipation rate correlations in the inertial 
range at equal times, sweeping effects do not matter. As we have shown previously |^] for 
the Burgers equation, for correlations involving velocities at different times sweeping effects 
can dominate the dynamic behavior: the structure functions, for example, S2 satisfy a wave 
equation with a characteristic velocity of the order of the rms fluctuations of velocity. 

While we have not proved that ^\ and \i2 are the only dominant exponents, especially in 
the case of the Navier-Stokes equation, our equations are exact, and the identifications we have 
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made in the 1-dimensional Burgers equation are transparent. The mathematical similarity of 
the equations derived from the Navier-Stokes equation to those in the Burgers case gives one 
confidence that one has indeed made progress in the derivation of the intermittency exponents 
in the Navier-Stokes case. 
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A Derivation of Eqns. (8) and (9) 



Starting from the Burgers equation for u\ = u{x\,t\) and multiply by u\ to obtain 

u\dui/dt\ = vu\d 2 u\jdx\ — \-du\jdx\ + u\T)\ . (20) 

The first term on the right-hand side can be rewritten using elementary manipulations as 

02 ip, 2 d(uidui/dxi) 2 

vu\d ui oxi = v\ — (du\ oxi) \ 

ox\ 

\&u\ldx\ - £1 . (21) 



Note that this is a key step since it allows us to disentangle the terms that are finite in the 
v — > limit, due to the singular behavior of the derivative terms, from those that vanish. 
Using this result and multiplying Eqn. (|2^) by e 2 yields 

2^r = " 12 " + e2Um (22) 

where we have used the fact that €2 does not depend on the coordinates x\ and t\. We 
compute next the average of the above equation with respect to the noise ensemble described 
by Eqn. (§). We assume that we have a spatially homogeneous and temporally steady state 
which implies that 

d < ... > /dt = and d < ... > /dx = . (23) 

Using the definitions x\ = x + r/2 and X2 = x — r/2, rudimentary calculus yields 

11 11 

d/dxi = -d/dx + d/dr = -d x + d T and d/dx 2 = -d/dx — d/dr = -d x — d r (24) 

and similar results for the time variable t\ and £2- Averaging Eqn. (^) yields 
^d T < e 2 u\ > = V -d 2 r < e 2 u\ > - < e x e 2 > - ^d r < e 2 u\ > + < e 2 «i?7i > . (25) 

Writing a similar equation with the subscripts 1 and 2 interchanged and adding yields the 
symmetric form displayed in Eqn. (||). 

The equation for the dissipation rate correlations in the Navier-Stokes case Eqn.(^) can 
be obtained analogously taking care to keep track of the Cartesian coordinate subscripts and 
the pressure term. One finds instead of Eqn. (|25"D, changing to primed variables as is done 
in the main text, the following equation 

-d T < e'u ■ u > = — V r 2 < e'u ■ u > +uV ri V rj < e'uiUj > 
— < e'e > —-d rj < e'u ■ uuj > 



2 



+ < e'u-fj> -d n < e'pu-i > (26) 
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The equations for the Burgers and the Navier-Stokes case correspond to each other term by 
term ( if one makes allowance for the fact that the one deals with one dimensional the other 
with three dimensional quantities) as is evident from the comparison of equations ( f25| ) and 
(j26[). The two additional terms in the Navier-Stokes case are due to the presence of pressure 
and the use of the full three dimensional definition of e, rather than its expression along one 
direction which is the only quantity usually measured Q. 



B Derivation of Eqn. (16) 

The general derivation of the equation satisfied by S q (r, r) is most easily carried out using 
generating functions. We use the usual notation u\ = u(x%,ti) and 112 = u(x2,t2), x\ = 
x + r/2, etc. Consider the generating function expa(iii — U2) and compute its derivative 
with respect to r. Starting from the stochastic Burgers equation for u\ and 112 and using 
d T = (1/2) (dt 1 — dt 2 ), we obtain 

= -(d tlUl + d t2 U2)e a ^- u ^ 

= ~e aiui - U2 \vdfui - uidwt + m + vdiuz - u 2 d 2 U2 + m) ■ (27) 
The following identities that are easily verified are useful: 

e a( Ul -u 2 ) Ui q iUi = 1 0T e a(ui-«a)( aUl _ ^ 1 

a 1 

e a( Ul -« 2 ) d 2 Ui = hdle aiui - U2) -a 2 e aiui - U2) (d lUl ) 2 }. (28) 

a 

We use these and the corresponding equations for derivatives with respect to x 2 and take 
average values over the Gaussian noise ensemble. We use Eqn. ( p4[) and assuming a 
spatially homogeneous, temporally steady-state which implies Eqn. ( p3[) we can evaluate 
< de a ^ Ul - U2 ^/dT > and obtain 

(eOCui-ua)) = _ 1^ (_^( e aK-« 2 ) (ni+U2)) + (e a( Ul - U2 ) (ei _ e2)) 

+ ^(e a{ " 2) (vi-V2)) ■ (29) 



The coefficient of a yields Eqn. (|iq) used in Sec. 3.3. 



C Galilean invar iance of Eqn. (9) 

In this Appendix we show the Galilean invariance of the righthand side of Eqn. (||) explicitly. 
A similar derivation can be carried out in the Navier-Stokes case. In the one-dimensional 
problem the Galilean transformation to a frame moving with a relative velocity V is accom- 
plished by letting 

x — > x + Vt and t —*■ t 

and 

d/dx — > d/dx and dt — > d/dt — Vd/dx. 
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Analogous results apply for the relative coordinates r = x\ — x 2 and r = t\ — t 2 . Making 
these substitutions in the righthand side of Eqn. (0) we find that terms proportional to V 3 
and V 2 vanish as can be easily checked. The term of order V is given by 

Id Id 1 v d 2 

2Q^< e i«2 - e 2 U! > +-— < ei«2 - e 2 uj > +- < e im + e 2m > + ^ ^ < e i n 2 + e 2 ui > . 

We now demonstrate an identity that shows that this term vanishes. Multiplying Burgers 
equation for u\ by e 2 we have 

d[e 2 ui]/dti = vd 2 [e 2 ui/}dx\ - ^d[e 2 u\]/dxi + uxr/i . (30) 

Adding the analogous equation with the subscripts 1 and 2 interchanged to the above equation 
and using Eqns. yields the required identity immediately. 



D Derivation of Eqn. (19) 



We derive the equation satisfied by < (Su ■ Su) 2 > where Su = u(x,t) — u(x',t') = u — u'. 
We use as usual 

x = R H — f and x' = R f 

2 2 



and 



1 l 

t = T +-T and t' = T r 

2 2 



leading to the identities 

0/dxi = ^d/dRi + d/dr; and d/dt = ^d/dT + d/dr 

and the corresponding equations for the primed variables. Using these and the Navier-Stokes 
equations for U{ and u\ and adding them yields the basic equation 



+ (3i) 

where 5p = p(z, t) — p(x', i'). Multiplying this equation by Su ■ Su Sui, summing over i and 
taking averages over the spatially homogeneous, temporally steady state of fully developed 
turbulence leads to the equation we wish to establish after straightforward manipulations. 
Terms that are proportional to v need careful consideration. It is most convenient to extract 
the terms that are finite as v — > by using Leibniz's rule: 

< Su ■ Su 5uiS7 2 5ui > = dj < Su ■ Su SuidjSui > 

— < Su ■ Su dj5uidj5ui > —2 < SuidjSuiSukdjSu^ > (32) 

and similar equation for the primed variable. The following identity that can be derived 
easily from incompressibility and Navier-Stokes equation 

- t £ < 33 > 
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and the obvious result that follows from it 

e = -^(diUj + djUi) 2 
are useful in rewriting the viscous terms. 



= u(diUj) 2 — vV 2 p 
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